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One of the most interesting properties of QCD is the existence of the deconfining 
transition at some temperature Tc, above which hadrons no longer exist and the 
strongly interacting matter is described in terms of quarks and gluons. Therefore 
an important task for lattice QCD is to study the nature of this transition and to 
O , determine the precise value of the transition temperature. In the next section I will 

discuss this problem more in detail. 

Although light hadrons cannot exist above the deconfinement temperature the 
situation can be different for heavy quarkonia. Because of the heavy quark mass 
k>( \ quarkonia binding can be understood in terms of the static potential. General 

^ ■ considerations suggest that quarkonia could melt at temperatures above the de- 

\ confinement temperature as a result of modification of inter-quark forces (color 

screening). It has been conjectured by that melting of different quarkonia states 
due to color screening can signal Quark Gluon Plasma formation in heavy ion col- 
lisions [Q. Many studies of quarkonia dissolution rely heavily on potential models 
[ 12 01 EI El El El ■ However it is very unclear if such models are valid at finite 
temperature [El. 

The problem of quarkonium dissolution can be studied more rigorously in terms 
of meson (quarkonium) spectral functions. Lattice calculation of charmonium spec- 
tral functions appeared recently and suggested, contrary to potential models, that 
J/tp and rjc survive at temperatures as high as 1.6Tc [0 E3 ^1- It has been 
also found that Xc melts at temperature of about l.lTc [^J 1121 113j . There are also 
preliminary calculations of the bottomonium spectral functions [ I14L I15| . In section 
3 I will discuss charmonium correlators and spectral functions. Finally conclusions 
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will be given in the last section. 

2. The finite temperature transition in QCD 

We would like to know what is the nature of the transition to the new state of 
mater and what is the temperature where it happens ^. In the case of QCD with- 
out dynamical quarks, i.e. SU(3) gauge theory these questions have been answered. 
It is well established that the phase transition is 1st order [EI- Using standard 
and improved actions the corresponding transition temperature was estimated to 
be Tc/v^ = 0.632(2) [[TTI (cr is the string tension). The situation for QCD with 
dynamical quarks is much more difficult. Not only because the inclusion of dynam- 
ical quarks increases the computational costs by at least two orders of magnitude 
but also because the transition is very sensitive to the quark masses. Conventional 
lattice fermion formulations break global symmetries of continuum QCD (e.g. stag- 
gered fermion violate the flavor symmetry) which also introduces additional com- 
plications. Current lattice calculations suggest that transition in QCD for physical 
quark masses is not a true phase transition but a crossover ^1 [201 1^ EH ■ The 
transition appears to be first order only for very small quark masses corresponding 
to pion mass of about 67MeV for three degenerate flavors [El- At finite temper- 
ature there is only one transition in QCD. The deconfinement transition identified 
with rapid increase in the degrees of freedom or increase of the Polyakov loop [ I23| 
coincides with the chiral transition [I18|. 

Recent lattice results for the transition temperature Tc from Wilson fermions 
[ 1^ 1^ . improved [El 1^ 1121 ^'^d unimproved staggered fermions [EDI with 2 
and 2-1-1 flavors of dynamical quarks are summarized in Fig. ^ The errors shown 
in Fig. n are only statistical with the exception of the data point from the MILC 
collaboration, where the large error partly comes from the continuum extrapolation 
and also includes systematic error in scale setting. One should note that calculations 
with Wilson fermions are done at very large quark masses (of the order of the strange 
quark mass or larger) and therefore the uncertainty in the value of Tc is dominated 
by the uncertainty of the extrapolation in the quark mass [jJSl and shown in Fig. ^ 
as an horizontal arrow. In the figure I also show the chemical freeze-out temperature 
at the highest RHIC energy [OH. Since the "critical" energy density ec = e(T'c) (i.e. 
the energy density at the transition) scales as the error in is the dominant 
source of error in [ I27j . 

3. Charmonia correlators and spectral functions 

In lattice QCD we calculate correlators of point meson operators of the form 

Jnit, x) = q{t, x)THq{t, x), (1) 

will talk here about the QCD finite temperature transition irrespective whether it is a true 
phase transition or a crossover and Tc will always refer to the corresponding temperature. 
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Fig. 1. Summary of lattiee results on the transition temperature Tc taken from 
Refs. [[H EHI EI El muni as well as unpublished result from S. Katz presented 
on Quark Matter 2005. The horizontal line and the band is the chemical freezout 
temperature at the highest RHIC energy [0^. 



where Th = 1, 75, 757/i, 7^71^ and fixes the quantum number of the channel to 
scalar, pseudo-scalar, vector, axial-vector and tensor channels correspondingly. The 
relation of these quantum number channels to different meson states is given in Tab. 
n Most dynamic properties of a finite temperature system are incorporated in the 
spectral function. The spectral function a{po,p,T) for a given mesonic channel H 
in a system at temperature T can be defined through the Fourier transform of the 
real time two-point fmictions D'' [xq, x, T) and (xq, x, T) or, equivalently, as the 
imaginary part of the Fourier transformed retarded correlation function [ '28'. The 
Euclidean time correlator calculated on the lattice 

Gh{t,P, T)^ I d\e'P-^TrJH{T, f) Jh(0, 0)) (2) 

is an analytic continuation of the real time correlator Gh{t,p^T) — (~iT,p,T). 
Using this fact and the Kubo-Martin-Schwinger (KMS) condition [ 28j for the corre- 
lators D^{xo,x,T) = D'^{xo + i/T,x,T), one can relate the Euchdean propagator 
Gh (t, p, T) to the spectral function through the integral representation 



G(T,p,r) - / du;aiu;,p,T)K{Lu,T,T) 







K(, , X T\ cosh(^(T - l/2r)) 

^("'"'^^ = sinh(c./2r) ■ 



To reconstruct the spectral function from the lattice correlator G{t, T) this integral 
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Table 1. Meson states in different channels 



representation should be inverted. Since the number of data points is less than the 
number of degrees of freedom (which is O(IOO) for reasonable discretization of the 
integral ) spectral functions can be reconstructed only using the Maximum Entropy 
Method (MEM) [ EH] • In this method one looks for a spectral function which 
maximizes the conditional probability P[(7\DH] of having the spectral function a 
given the data D and some prior knowledge H which for positive definite spectral 
function can be written as 

P[a\DH]=cM-lx^+c^S), (4) 

where 5* = J dLj[a(Lj) — 'm{LLj) — a{LL>)hi{a{LL>)/'m{uj))] is the Shannon - Janes entropy. 
The real function m{u!) is called the default model and parametrizes all additional 
prior knowledge about the spectral functions, such as the asymptotic behavior at 
high energy [I29|. In order to have sufficient number of data points either very fine 
isotropic lattices [E] ^1 El or anisotropic lattices [0 El El have been used. 

The spectral function for pseudo-scalar charmonium spectral functions calcu- 
lated on anisotropic lattice [El is shown in Fig. [2 The first peak in the spectral 
function corresponds to ric{lS) state. The position of the peak and the correspond- 
ing amplitude (i.e. the area under the peak) are in good agreement with the results 
of simple exponential fit. The second peak in the spectral function is most likely 
the combination of several excited states as its position and amplitude is higher 
than what one would expect for pure 2S state. The spectral function becomes sen- 
sitive to the effects of the finite lattice spacings for w > 5GeV. In this uj region the 
spectral functions becomes also sensitive to the choice of the default model. Also 
shown in Fig. [3 is the spectral function in the scalar channel from Ref. [El- The 
1st peak corresponds to Xcoi^P) state. The correlator is more noisy in the scalar 
channel than in the pseudo-scalar one. As the results the Xcoi^P) peak is less pro- 
nounced and has larger statistical errors. The peak position and the area under 
the peak is consistent with the simple exponential fit. As in the pseudo-scalar case 
individual excited states are not resolved and the spectral function depends on the 
lattice spacing and default model for lu > 5GeV. Similar results have been found 
for the vector and axial-vector channels which correspond to J/^p and Xci states 
respectively. 

We would like to know what happens to different charmonia states at tem- 
peratures above the deconfinement temperature Tc. With increasing temperature 
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Fig. 2. Charmonium spectral function in the pseudo-scalar channel (left) and the 
scalar channel (right) at different lattice spacings and zero temperature from Ref. 
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Fig. 3. The ratio G{t,T) /GreconiT,T) of charmonium for pseudoscalar channel 
at at a^^ = 14.11GeV (left) and scalar channel at at a^^ = 8.18GeV (right) at 
different temperatures [El- 



it becomes more and more difficult to reconstruct the spectral functions as both 
the number of available data points as well as the physical extent of the time di- 
rection (which is 1/T) decreases. Therefore it is useful to study the temperature 
dependence of charmonia correlators first. From Eq. Q it is clear that the temper- 
ature dependence of charmonia correlators come from two sources: the temperature 
dependence of the spectral function and the temperature dependence of the inte- 
gration kernel K{t,u!,T). To separate out the trivial temperature dependence due 
to the integration kernel, following Ref. for each temperature we calculate the 
so-called reconstructed correlator Grecon{T, T) = diua{u!, T = 0)K{t, uj, T). Now 
if we assume that there is no temperature dependence in the spectral function, then 
the ratio of the original and the reconstructed correlator should be close to one, 
G{T,T)/Grecon{T,T) ~ 1. This way we can identify large changes in the spectral 
function. This gives reliable information about the fate of charmonia states above 
deconfinement. In Fig. Owe show this ratio for pseudo-scalar and scalar channels 
correspondingly calculated on anisotropic lattice From the figures one can see 

that the pseudo-scalar correlators shows only very small changes till l.STc indicating 
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that the rjc state survives till this temperature with little modification of its proper- 
ties. On the other hand the scalar correlator shows large changes already at l.lGTc 
suggesting strong modification or dissolution of the Xco state at this temperature. 

More detailed information on different charmonia states at finite temperature 
can be obtained by calculating spectral functions using MEM. The results of these 
calculations are show in Figs. 0] Because at high temperature the temporal ex- 
tent and the number of data points where the correlators are calculated become 
smaller the spectral functions reconstructed using MEM are less reliable. To take 
into account possible systematic effects when studying the temperature modifica- 
tions of the spectral functions we compare the finite temperature spectral functions 
against the zero temperature spectral functions obtained from the correlator using 
the same time interval and number of data points as available at finite temperature. 
We see that spectral function in the pseudo-scalar channel show no temperature de- 
pendence within the statistical errors. This is in accord with the analysis of the 
correlation functions. Also the spectral functions show very little dependence on 
the default model. Similar conclusion has been made in Ref. [ 1111 IT^ where cor- 
relators and spectral functions have been calculated on very fine isotropic lattices 
as well as in Ref. [ lO where anisotropic lattice have been used. The pseudo-scalar 
spectral function was found to be temperature independent also in Ref. [ 9 where 
correlators of extended meson operators have been studied on anisotropic lattices. 
The study of the charmonium correlators with different spatial boundary condi- 
tions provides further evidence for survival of the IS" charmonia states well above 
the deconfinement transition temperature [ 30 . 

The scalar spectral function shows large changes at l.lGTc which is consistent 
with correlator-based analysis. Also the default model dependence of the scalar 
correlator is large above the deconfinement transition (c.f. Fig. 3, right). This 
means that the XcO i'^Po) dissolves at this temperature. Similar results for the 
scalar spectral function have been reported in The results for the axial 

correlators and spectral functions are similar to scalar ones as expected. 

The vector correlator, however, has temperature dependence different from that 
of the pseudo-scalar channel [ 31 . This is due to the fact that the vector current is 
conserved and there is a contribution to the spectral functions at very small energy 
w ~ corresponding to heavy quark transport [ ,32. .33,. The transport peak in the 
spectral functions can be written as [ \'62\ 



T 1 UJ7] 

^^oM-XsiT)---^,, (5) 



where r] = T/AI/ D with D being the heavy quark diffusion constant. Furthermore 
Xs{T) is the charm or beauty susceptibility and M is the heavy quark mass. To the 
first approximation the transport contribution to the spectral function gives rise 
to a positive constant contribution to the correlator Giow{t) — Xs{T)T'^/M [I32|. 
resulting in the enhancement of the finite temperature correlators relative to the 
zero temperature ones, in agreement with the lattice data presented in [ 1311 I34j . 
Finite value of the diffusion constant D will give rise to some curvature in Giow (t) . 
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Fig. 4. Charmonia spectral function in the psuedoscalar channel at a^^ = 
14.11GeV (left) and the scalar channel (right) at a^^ = 8.18GeV for zero and 
finite temperature [I13|. For finite temperature scalar channel two different default 
models are shown. 

The smaller the value of D is, the larger is the curvature in Giow{t). Thus extracting 
G/oto(T) from lattice data and estimating its curvature can give an estimate for D. 
This, however, requires very precise lattice data which are not yet available [I32|. 

4. Conclusions 

In this contribution I discussed the status of finite temperature lattice QCD cal- 
culations. The transition in QCD appears to be a crossover and not a true phase 
transition. The value of the transition temperature, however, is not well know. 
Current lattice calculations give estimates between 167 and 188 MeV with large 
statistical and systematic errors. Recently charmonium correlators and spectral 
functions have been calculated in lattice QCD and show that IS charmonia states 
survive above the transition temperatures till 1.6Tc while the IP state dissolve at 
T ~ l.lTc. Furthermore, the properties of the IS" charmonia states are not modified 
siginificantly above Tc- On the other hand potential model predict strong modifica- 
tion of charmonia properties and therefore are not consistent with the lattice data 
on the correlators presented above [I33|. 
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